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Abstract. We study the dynamics of discrete-time regulatory networks on random 
digraphs. For this we define ensembles of deterministic orbits of random regulatory 
networks, and introduce some statistical indicators related to the long-term dynamics of 
the system. We prove that, in a random regulatory network, initial conditions converge 
almost surely to a periodic attractor. We study the subnetworks, which we call modules, 
where the periodic asymptotic oscillations are concentrated. We proof that those modules 
are dynamically equivalent to independent regulatory networks. 
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1. Introduction 

Numerous natural and artificial systems can be thougli as a collection of basic units in- 
teracting according to simple rules. Examples of this interacting systems are the ge- 
netic regulatory networks, composed of interactions between DNA, RNA, proteins, and 
small molecules. In social or ecological networks, a similar regulatory dynamics may 
also be considered. The traditional way to model these systems is by using coupled 
differential equations, and more particularly systems of piecewise affine differential equa- 
tions (see [6, 9, 17]). Finite state models, better known as logical networks, are also used 
(see [11, 13, 19]). Within these modeling strategies, the interacting units have a regular 
behavior when taken separately, but are capable to generate global complex dynamics 
when arranged in a complex interaction architecture. In all the models considered so far, 
each interacting unit regulates some other units in the collection by enhancing or repress- 
ing their activity. It is possible then to define an underlying network with interacting 
units as vertices, and their interactions as arrows connecting those vertices. The theo- 
retical problem we face here is to understand the relation between the structure of the 
underlying network, and the possible dynamical behaviors of the system. We will do this 
in the context of a particular class of models first introduced in [22] , and further studied 
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in [5] and [16]. In these models, the level of activity of each unit is codified by a positive 
real number. The system evolves synchronously at discrete time steps, each unit following 
an affine contraction dictated by the activity level and interaction mode of its neighboring 
units. The contraction coefficient of those transformations determines the degradation 
rate at which, in absence of interactions, the activity of a given unit vanishes. In the 
framework of this modeling we have proved general results concerning the constrains im- 
posed by the structure of the underlying network, over the possible asymptotic behaviors 
of a fixed system [16]. In the present paper, following [22], we will focus on the asymptotic 
dynamics of regulatory systems whose interactions are chosen at random at the beginning 
of the evolution. Within this approach, individual orbits are elements of a sample space, 
and the statistical indicators we will study become orbit dependent random variables. The 
probability measures we used are built from a fixed probability distribution over the set 
of possible underlying networks. Then, given a fixed underlying network, we associate a 
sign to each one of its arrows, depending on whether the interaction they represent are 
activations or inhibitions. Positive and negative signs are randomly chosen, keeping a 
fixed proportion of negative arrows inside a given statistical ensemble of systems. In this 
way it is possible to study certain characteristics of the asymptotic dynamics, as function 
of the degradation rate and the proportion of inhibitory interactions. 

Our first result states that in a random regulatory network, initial conditions converge 
almost surely to a periodic attractor. This result points to the conclusion that in regulatory 
dynamics, the origin of the complexity is the coexistence of multiple dynamically simple 
attractors. We prove that this is the case in a full measure set in the parameter space. 

According to our preliminary numerical explorations, the long-term oscillations of the 
system concentrate on subnetwork whose structure depends on the parameters of the 
statistical ensemble of regulatory networks. Our second result states that the dynamics one 
can observe when restricted to the oscillatory subnetwork, is equivalent to the dynamics 
supported by the subnetwork considered as an isolated system. This result allows us 
to introduce the concept of modularity. If we call module any observable oscillatory 
subnetwork, then, according to our result, the dynamics of a small network is preserved 
when it appears as a module in a larger network. We interpret this as the emergence of 
modularity. This result allows us to predict admissible asymptotic behaviors in regulatory 
networks admitting disconnected oscillatory subnetworks. It is worth mentioning that 
this kind of modularity was already studied in the context of Boolean networks [3] and 
more recently, in continuous-time regulatory networks [10]. Our approach allows a formal 
approach to the problems addressed in those works. 

The paper is organized as follows. In the next section we will introduce the objects under 
study, then in Section 3 we will state the results, and present some of the proofs. After 
reviewing two examples, which we do in Section 4, we will give the proofs of the two 
more technical results in Section 5. The paper ends with a section of final comments and 
conclusions. 

This work was supported by CONACyT through the grant SEP-2003-C02-42765, and by 
the cooperation ECOS-CONACyT-ANUIES M04-M01. E. U. thanks Bastien Fernandez 
for his suggestions and comments. 
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2. Preliminaries 

2.1. Regulatory Networks as Dynamical Systems. 

The interaction architecture of the regulatory network is encoded in a directed graph 
Q = {V, A), where vertices V represent interacting units, and the arrows A C V xV denote 
interaction between them. To each interaction (n, v) A we associate a threshold T(^u,v) S 
[0, 1], and a sign ai^u,v) ^ {"Ij 1} that is chosen according to whether this interaction is an 
inhibition or an activation. We quantify the activity of each unit v (zV with real number 
x„ G [0, 1]. Thus, the activation state of the network at a given time t is determined by 
the vector x* € [0,1]^. The influence of a unit u over a target unit v turns on or off, 
depending on its sign, when the value of x„ trespasses the threshold T(^u,v) ■ The evolution 
of the network is generated by the iteration of the map Fg^cr,T,a '■ [0, 1]^ — > [0, 1]^ such 
that 

(1) x*+i = Fe,,,T,a(x*) := ax* + (1 - a)Z5g,,,T(x*), 

where the contraction rate a G [0, 1) determines the speed of degradation of the activity 
of the units in absence of interaction, and the interaction term Z)g^o-,T : [0, 1]^ [0, 1]^ is 
the piecewise constant function defined by 

(2) £>g,^,T(x)„ := H{a{^^^^) (x„ - T(^u,v))) ■ 

Here : R — > M is the Heaviside function, and Id(f) := ^ V : {u,v) S A} stands 
for the input degree of the vertex v. We will be referring to the discontinuity set of the 
transformation Dg^^^T, which is 

(3) At := {x e [0, 1]^ : x^ = Tu^v for some(u, e A^ . 

For each a G [0, 1), the transformation Fg^„^T,a is a piecewise affine contraction with 
discontinuity set A^. 

From now on, by a discrete-time regulatory network we will mean a discrete-time dy- 
namical system ([0, 1]^, Fg^o-,T,a)) with phase space [0,1]^, and evolution generated by 
the piecewise affine contraction Fg^a-,T,a '■ [0)1]^ — *■ [0;1]^ defined in Equation (1). The 
discrete-time regulatory networks studied here have interactions of equal strength, and 
they act additively on each target unit. More general discrete-time regulatory networks 
have been considered in [5, 16]. 

2.2. Statistical Ensembles. 

We build our statistical ensembles as follows. We fix the value of the contraction rate 
o G [0, 1) and the set V representing the interacting units. Then we consider the set of all 
possible piecewise transformations, 

(4) ^a,y := {Fg,<x,T,a : g:={V,A), aG{-l,l}^, Tg[0,1]^}. 

The individual elements of our statistical ensembles are couples (Fg^o-,T,a5 x), with Fg^fj,T,a £ 
Tay and x G [0,1]^. A couple (Fg^o-.T.a, x) determines a deterministic orbit {x* := 
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We supply the sample space Tay ^ [0, 1] with a probability measure as follows. First we 
fix a probability distribution ¥g over the set C5y := {(y, j4) : A C V x V} of all directed 
graphs with vertex set V. Then, for i] E [0, 1] and Q = {V,A), we choose sign —1 with 
probability r] and +1 with probability 1 — 7], independently for each arrow in Q. The 
thresholds are independent and uniformly distributed random variables in [0, 1]^, as well 
as the initial conditions in [0, 1]^. In this way we obtain the probability measure Fa^n on 
Ta,v X [0,1]^ such that 



(5) Pa,„ {(Fg,,,T,a,x) : T£l,^£j} = FaiQ) X Pa,„(c7) X vol(T) X Y0\{J), 

for all rectangles I C [0, 1]^ and J C [0, 1]^. Here Pa,,, : {-1, 1}^ ^ [0, 1] is such that 



{u,v)£A ^ 

and vol stands for the Lebesgue measure on the corresponding euclidean spaces. 

2.3. Some graph theoretical notations and definitions. 

A path in ^ = (^4, V) is a sequence uq, ui, . . . ,Uk, with Ui £ V and (uj, ttj+i) G A for each 
i. The length of a path is the number of arrows it contains. A cycle is a path uo,ui, . . . ,Uk, 
where uq = u^- We say that the vertices u,v £ V are connected if there exists a cycle 
uo,ui, . . . , Mfc such that u,v £ {uq,ui, . . . , Uk}- The connected components of Q are the 
maximal subgraphs of Q such that all their vertices are connected. The distance between 
two vertices u, v £ V, which we denote dg{u,v), is the length of the shortest directed 
path whose end vertices are u and v. 

2.4. The oscillatory subnetwork and the asymptotic period. 

To each couple (Fg^o-.T.a, x) we associate a directed graph ^osc := (^osci^osc) ^ G, the 
oscillatory subnetwork, defined by 

^osc := {{u,v) £ A : -ff(x^ — T^^v) does not converge}, 
V'osc := {v£V: 3 u £V such that {{u,v),{v,u)} n Aosc 7^ 

By definition, the oscillatory subnetwork is spanned by all the arrows whose activation 
state H{au,vi^u — Tu,v)) changes infinitely often. This subnetwork is the equivalent, 
in discrete-time regulatory networks, to the dynamical islands introduced in [10] for 
continuous-time regulatory networks. They are also related to the clusters of relevant 
elements in Boolean networks, as they were defined in [3]. Unlike the differential equa- 
tions and the finite state modeling strategies of regulatory dynamics [14, 20, 21], our 
models admit oscillations in all network topologies, regardless the distribution of activa- 
tions and inhibitions through the graph. This is possible because of the discreteness of 
time, and continuity of the state variable. Therefore, oscillatory subnetworks may occur 
in any discrete-time regulatory network. 

The oscillatory subnetwork can be seen as a random variable 



(6) 




^osc : J'ay X [0, 1] 



V 



(Scv ■■= {{V, A) :V CV and AcV xV} 
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from which we can derive other statistical indicators, as for instance its size i^Vosc, the 
number nc(t/osc) of its connected components, and its the degree distribution 

#{v G Vosc : Idosc(u) = k} 

PGosAk) := — , 

if Vosc 

where Idosc('v) := {u G Ksc : {u,v) G Aqsc}- 

Another statistical indicator we will consider is the asymptotic period of a given orbit. 
Denote by PerT-(Fg o- ^- a) the set of all o- T- Q-periodic points of minimal period r, i. e., 

(7) Per,(Fg,,,T,a) = {x € [0, 1]^ : FgV,T,a(x) = x and F^g^^^^M / x if t < r} . 

We will say that the asymptotic Fg^o-^y^a^period of x G [0, 1]^ is r if there exists y G 
Y'evT-{Fg^a^T,a): such that 

lim |i^kr,a(x)-4.,r,a(y)l=0. 

t— ♦OO ' ' »j 1 1 

We will denote this by P(Fg^o-,T,ai x) = r. 



3. Results 



3.1. The asymptotic period. 

To each finite set V , and a G [0, 1), we associate the sample space Tay x [0, 1]^, with Tay 
as in Equation (4). We supply this sample space with the product sigma-algebra, taking 
for [0, 1]^ and [0, 1]"^, the corresponding Borel sigma-algebras. The next result concerns 
the asymptotic period (Fg;,^,T,a, x) i-^ P{Fg^„^T,a,^)- 

Theorem 1. Given a finite set V , and a G [0, 1), the function P : Tay x [0, 1]^ N 
which assigns to (Fg^o-,T,a, x) G Tay x [0,1]^ the value of the asymptotic Fg^cr,T,a-period 
of X, is a measurable function. Furthermore, for each rj G [0, 1], 

Pa,r, {(Fg,,,T,a,x) G J^ay X [0, 1]^ : P(Fg,,,r,a, x) < CX)} = 1. 

This expected but nontrivial result says that a random orbit will almost certainly approach 
a periodic orbit, otherwise said, an initial condition almost certainly converges to a periodic 
attractor. It is because of this result that we can restrict ourselves to the study of periodic 
orbits, even though the parameter set for which the orbits are infinite may be uncountable 
(see [5]). We could therefore redefine the oscillatory subnetwork by saying that, after a 
transitory, their arrows change periodically their activation state. The proof of this result 
and some related comments are left to Subsection 5.1. 
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3.2. Modularity. 

In this paragraph we wih consider probabihty measures obtained by conditioning on a 
subnetwork, from a given statistical ensemble of regulatory networks. Given ¥a^n on J-a,v ^ 
[0, 1]^, and a subnetwork Q := {V, A) G ©cy, we define the analogous probability measure 



g on Tg^Y X [0, 1]^, such that 



(8) Pa,,,c7 { (^g,^,T,a, y)-TeI,yeJ}.-^^ " otherwise, 

for each a G {—1,1}^, and all rectangles 2 C [0,1]"^ and J C [0,1]^. As before, vol 
denotes the corresponding Lebesgue measures. This measure can also be seen as the 
marginal of Pay obtained by projecting over V C V, then conditioned to have underlying 
network Q & '3y. 



Theorem 2. Fix a G [0, 1), r] G [0, 1] , and a digraph Q = {V, A) G ©cy with V CV. If the 
digraph distribution Fq is such that FciQ) > for all Q G ©y, and i/Pa,7?(^osc = Q) > 0, 
then we can associated to Q: 

a) a digraph extension ^ext := iV,A) G ©y, 

b) rectangles I := I x I' C [0, 1]^ x [0, and J := J x J' C [0, 1]^ x [0, , 

c) and for each a G {—1, 1}"^; o-n extension cjcxt G 1}^ such that ciextU ~ ^■ 
These associated objects satisfy 

^osc(Fg-,,„^,,„T,a>x) <ZG MT(^T,i^(^J, anda(^ {-1, 1}^. 
Furthermore, Qcxt, T, and J , define surjective transformations 

$A : ^ X 2:' ^ [0, 1]^, such that <^a{T x T') = V^f + with diagonal and 
constant, 

e) and similarly, ■ >J x >7' ^ [0, 1]^ such that <l>\/(x x x') = PyX + Cy, with T>y 
diagonal and Cy constant. 
These transformations satisfy 

for all T and x G ^7. 

The proof of this result follows from a construction, and it is postponed to Section 5. Let 
us from now expose its interpretation and some of its consequences. We can divide this 
Theorem into two parts, the first part concerns the stability of the oscillatory subnetworks. 
It says that if an oscillatory subnetwork Q := {V, A) has positive probability to occur, 
then to each signs matrix a G {—1, 1}"^ it corresponds an oscillatory subnetwork which 
is subgraph of Q. Furthermore, the set of initial conditions and thresholds for which a 
given subgraph of Q is the oscillatory subnetwork, is a rectangle with nonempty interior. 
Therefore, each of those subgraphs is stable under small changes in the thresholds T G 
[0,1]^, and in the initial condition x G [0,1]^. The size of the maximal perturbation 
depends on the position of (T, x) with respect to the borders of the above mentioned 
rectangle. 
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More interestingly, the second part of the theorem states that if an oscillatory subnetwork 
has positive probability to occur, then the dynamics one can observe when restricted to 
that subnetwork is equivalent to the dynamics supported by the subnetwork considered as 
an isolated system. The equivalence is achieved through the transformation and ^>y. It 
is because of this equivalence that we can talk about modularity. Indeed, if we call module 
any oscillatory subnetwork occurring with positive probability, then this theorem can be 
rephrased by saying that each orbit admissible in a module considered as an isolated 
system, can be achieved, up to a change of variables, as the restriction of an orbit in 
the original system. The network extension Gext, the values of the external thresholds 
T G [0, and the external initial conditions J' E [0, can be though as the 
analogous in our systems, to the functionality context defined in [18] for logical networks. 
Using their nomenclature, our theorem ensures that if an oscillatory subnetwork has a 
positive probability to occur, then there is positive measure set of contexts making this 
subnetwork functional. 

An interesting consequence of the previous theorem is the following. 

Corollary 1. Fix a G [0, 1), r] G [0, 1], and Q = Q1UQ2 ^ ©cy, with Qi := {Vi.Ai) and 
Q2 := (V2, A2) vertex disjoint and such that Vi U V2 C y. If the digraph distribution is 
such that Pg(^) > for all G G (Sy, and i/Pa,»?(^osc ^ ^) > 0, then 

Pa,^(P = t)-^C Y1 ^->vMP = ri) K,v,G2iP = ^2), 

lcm(ri ,T2)=T 

with C > a constant depending on Fa^^ and Q. 

The probabilities in the sum at right hand side of the inequality are analogous probability 
measures of the kind defined by Equation (8). They are obtained as marginal probability 
corresponding to projections over vertex the sets Vi and V2, then conditioned to have 
underlying network Qi and Q2 respectively. 

In [3] the authors study the relationship between the modular structure and the periods 
distribution in Boolean networks. The previous corollary establishes such a relation in 
the framework of discrete-time regulatory networks. Our result relates the distribution of 
the asymptotic period of a network Q, to the distribution of the asymptotic period of its 
oscillatory subnetworks considered as isolated systems. It establishes in particular that 
the observable asymptotic periods of Q are the least common multiples of the observable 
asymptotic periods of its oscillatory subnetworks. 

Proof. First notice that, for each Q G (5cy we have 

Let G = {V,A) = {Vi U V2,Ai U A2), with Gi := {Vi,Ai) and ^2 := (1^2, ^2_) vertex 
disjoint. Theorem 2 ensures the existence of the following objects associate to ^: a) the 
extension of G, a digraph ^cxt := (^)^) € &v such that G ^ Gcxt, b) two affine surjective 
transformations 

$A : X ^ [0, X [0, 1]^2 and ■ J ^ [0, 1]^^ x [0, lf\ 
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and c) for each a = a^^^ x a^"^^ G { — 1, 1}"^^ x { — 1, 1}^^, an extension cJcxt € { — 1, 1}"^ such 
that 0"|aiuA2 = ^- From the second part of the theorem it follows that 



(0 



for all t E N and (r,x) elxj. Here we have used $i/(T) = T^^) xT^^) g [0, x [0, 1]^^ 
and $i/(x) = y(i) x y(2) g [0, x {0,1]^^. For i = 1,2, and each a(^) € {-1,1}^^ and 
Tj G N, let us define 

Defining <^ : I x J ^ [0,1]^ x [0,1]^ such that $(T,x) := $a(T) x $y(x), it follows 
from (9) that 

lcm(ri ,T2)=T 

where vol denotes Lebesgue measure in [0, 1]"^ x [0, 1]^. Let us now use the decomposition 
1 = 1 xT C [0,1]^ X [0,1]^\^ and J = J X J' C [0,1]^ x [0,1]^\^. According to 
Theorem 2, we have 

$((r X T') X (x X x')) = V{f X x) + C 
where D is a diagonal linear bijection, and C is a constant. With this, 

F..(P = r|^osc = a) ^ ^^(^oxOyol(rx_^0 

X Yl (^gi,<^(l),ri) X V0I2 (^g,,<,(2),,J , 

lcm(Ti ,r2)=T 

where vol denotes the Lebesgue measure in [0, 1]"^^"^ x [0, 1]^^^, volj denotes the Lebesgue 
measure in [0, 1]"^* x [0, 1]^" for i = 1,2, and I'D] denotes the determinant of the transfor- 
mation D. From here, and taking into account the definition in Equation (8), we have 

IPG(acxt) vol(Z' X J') mm{rj, 1 - 7?)#(^\^) 



\n {P = tIGosc = G) > 



\V\ Pa,„(e?osc C Q) 



X > F^,.„ (£7^1)) X Pa2,,? (^^^^ 



j7=f7(l)Xf7(2) 

X 



^ voli [Pg^-aW.j^ X V0I2 (^g2,<j(2),, 

lcm(ri ,T2)=r 

IPG(gext) vol(3:^ X J') mm{7], 1 - 

lcm(Tl,T2) = T 
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and the result follows with C := PG(^cxt) x vol(2:' x J') x mm{r], 1 - r])*^^\^^ /\V\. □ 
3.3. Sign Symmetry. 

The next result illustrates how a structural constrain in the underlying random network 
manifests in the distribution of the statistical indicators. Here, the fact that the underlying 
random digraph does not admit cycles of odd length, implies a symmetry on the statistical 
indicators with respect to the change of sign in the interactions. Let us remind that for 
each rj G [0, 1] and a G [0, 1), the probability measure Pa,r? on Ta,v x [0; 1]^ which defines 
a statistical ensemble, is obtained from a fixed distribution Fq over the finite set of all 
directed graphs with vertices in V. We have the following. 

Proposition 1. // the distribution Fq over the set &v of all directed graphs with vertices 
in V is such that ¥g{Q G C5\/ : Q admits a cycle of odd length } = 0, then 

Pa,„(P = r) = Pa,i_^(P = r), VrGN, 

Pa,„(^osc = ^) = ^a,l-v {Gosc = G) V ^ G ©CV, 

for each rj G [0, 1] . 

Statistical ensembles of regulatory networks whose underlying digraphs are random trees 
(as defined in Subsection 4.1 below) clearly satisfy the hypothesis of Proposition 1. Ran- 
dom subnetworks of the square lattice also have this property. Therefore, for regulatory 
dynamics over those kind of digraphs, the statistical indicators we consider are left invari- 
ant under the symmetry 77 1-^ 1 — 77. 

In the proof of Proposition 1, we will need the following. 

Lemma 1. For each Q := (V, A) G <5v, and a G {—1, 1}^ we have, 

vol{(r,x) G [0,1]^ X [0,1]^ : {FI^tM): t G N} n At = 0} = 1. 

where At := {x G [0, 1]^ : Xu — Tu^i, for som,e(^u, f) G ^4} is the discontinuity set of the 
piecewise constant part of Fg^(j^T,a- 

This lemma directly follows from the arguments developed below, in paragraphs 5.1.3 
and 5.1.4, inside the proof of Theorem 1. 

Proof of Proposition 1. First notice that, for each Q = [A, V) G and a G { — 1, 1}"^ we 
have 

FA,i-,{-a) = n {^^^^^f^ + <^M{l-ri) 

= n 

{u,v)eA 

Let us suppose that all the cycles in Q have even length, and for each connected component 
of G ■= {A,V) ^ G = [A,V) G CJy, choose a pivot vertex u ^ V. With this define 



a 



{u,v] 



+ 1 



{u,v)'n 
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■■ [0, 1]^ ^ [0, 1]^ and ■ [0, 1]^ ^ [0, 1]^ such that 



_ I 1 - r(„^„) if dg(n, n) G 2N _ / 1 " x« if dg{u, u) € 2N 

• 1 T(„^^^ otherwise "^VK^ju ■ <| otherwise, 



for each u,v £ V. Here dg denotes vertex distance induced by the digraph G, as defined in 
Subsection 2.3. Both and are affine isometrics, therefore they preserve the Lebesgue 
measure m [0,1]^ and [0,1]^ respectively Udg{v, u) is even, then \I'y(x),u = 1 — x^, and 
for each u such that {u,v) € A, ^'\/(x)„ = x„ and ^a{T){^u,v) = '^{u,v)- I^i this case we 
have 

„^^^^(T),a(^\/(x))^ = a(l - x^) + jrp: ^ H (-cr(„,„) (x„ - T(„,„))) 

= «(i - x^-) + fdf^ (x„-r(„,„)))) 

for each x ^ A^- On the other hand, if dg{v^ I'pivot) is odd, we have (x)^, = x^,, and for 
each u such that {u,v) G yl, '^y{^)u = 1 — x^ and ^'a(T)(„^^) = 1 — T(^u,v)- In this case 

-a,vI'^{T),a(^v(x))t, = ax^ + ^^T^ J2 ^ {-^iu,v) ((1 - ^u) - (l - ^u,v)))) 

ueV:{u,v)eA 
^ ' ueV:iu,v)eA 

= -Fg,(7,T,a(x)t, = ^y(Fg,„,T,a(x))^, 

for each x At- Taking into account Lemma 1, Fg ^(^'v(x)) = ^(x)) 

for all t € N, for almost all (x, T) G [0, 1]^ x [0, 1]^. Since $y is an isometry, we immedi- 
ately have 

V0l{(r,x) G [0,1]^ X [0,1]^ : P(Fg,_,,q,^(r),a,*l/(x)) =P(Fe,,,T,a,x)} = 1, 

for all G := (V,A) e &v with no odd cycles, and every a G {—1, l}"^- Since both ^a and 
preserve the respective Lebesgue measure, and Pg('S admits an odd cycle) = 0, we 
obtain 

xvol{(r,x) : P(Fg,-.,*^(T),a,^y(x)) =t} 
= Y^^GiQ) Yl ff^A,i-r,(-fT) vol{(r,x) : P(Fg,_,,T,a,x) =r} 

= Pa,i_^(P = r), 

for all r G N, and this concludes the proof of the first claim in the proposition. 
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For the second claim, if Fg ^-y ^(x) ^ At, then 

H (a(„,,)(F*^,^r^„(x)„ - r(„,,))) = H {-au,v {^v{F^,a,T,ai''))u - ^a(T)(„,„))) 
Therefore, by Lemma 1, 

V0l{(r,x) G [0,1]^ X [0,1]^ : aosc(i^e,-<x,.I>^{T),a,^y(x)) = eosc(i^e,a,T,a,x)} = 1, 

for all := {V,A) G &v all of whose cycles have even length, and every a € {—1,1}"^. 
From here, taking into account that Pg(^ admits an odd cycle) = 0, and the fact that 
that and preserve the Lebesgue measure, we obtain 

Q&'&V o-G{-l,l}^ 

XV0l{(r,x) : aosc(i^e,-a,*^{T),a,^y(x)) 
= Y.^g{G) lPA,l-r,(-f^) V0l{(r,x) : gosc(i^e,-.,T,a,x) 

(tg{-i,i}-4 

= Pa,l-;,(6^osc = 

for all Q € (Scy, and the proof is completed. 

□ 

4. Examples 

In this section we present, as a matter of illustration, two families of random digraphs 
which we have explored numerically. A detailed numerical study, which requires massive 
calculations, is out of the purpose of the present work, and is left for a future research. 
Instead, in this section we make some general observations suggested by a preliminary 
numerical exploration of these examples. 

4.1. Two Families of Random digraphs. 

We consider two kind of distributions Pg on the set 0y := {(y, A) : A dV y~ V} of all 
directed graph with vertex set V . On one hand we have a directed version of the classical 
Erdos-Renyi ensemble, which we define as follows. 

Fix a vertex set V and p S (0, 1), and consider the probability distribution Pp : ^ [0, 1] 
such that Pp(^) = p*"^. According to this, a directed graph contains the arrow {u,v) 
with probability p, and the inclusion of different vertices are independent and identically 
distributed random variables. The typical Erdos-Renyi graph is statistically homogeneous, 
thus suited for a mean field treatment. Most of the rigorous results concerning random 
graphs refer to these kind of models (see [4, 8] and references therein). 

The other family of random digraphs we will refer to derives from the famous Barabasi- 
Albert model of scale-free random graphs, whose popularity relies on the ubiquity of the 
scale-free property [2, 15]. A dynamical construction of scale-free graphs was proposed 
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by Barabasi and Albert in [1] (see [8] for a more rigorous presentation). Their model 
incorporates two key ingredients: continuous growth and preferential attachment. We 
implemented their construction as follows. For n = 0, we take Qq = Kmo, the complete 
simple undirected graph in mo vertices. Then, for each n ^ 1, a new vertex Vn+i is added 
to the graph Qn ■= (EmVn)- This new vertex form new edges with randomly chosen 
vertices in Vn- The probability for v £ Vn to he chosen is proportional to its degree in Qq, 
i. e, 

{vn+i,v\ £ En With probability Pn{v) := —= . 

At the {n + l)-th step we obtain a graph Qn+i with an increased vertex set Vn+i ■= 
Vn U and an enlarged edge set En+i- This random iteration continues until a 

predetermined number N of vertices is obtained. In this scheme we are able to add more 
than one edge at each iteration. If on the contrary, we allow only one new edge at each 
time step, and we start with mo = 2 vertices choosing at the 1-th step one of the two 
preexisting vertices to form a new edge with with probability 1/2, then the resulting graph 
would be a random scale-free tree [8]. In any case, the final graph GN-mo is turned into a 
directed graph Q = iV, A) by randomly assigning a direction to each edge {n, v} £ EN^nio ■ 
In our case we choose one of the two possible directions with probability 1/4, and give 
probability 1/2 to the choice of both directions. This procedure generates a probability 
distribution Pba in ■= {(^,^) -.AcVxV}. 

4.2. Erdos Renyi. 

We generated random digraphs of the Erdos-Renyi kind with 100 vertices, and probability 
of connection p £ {0.2,0.4,0.6,0.8}. For p fixed and each ry £ {0, 0.1, 0.2, . . . , 1}, we 
build 20 graphs, to which we randomly and uniformly assign thresholds, and signs with 
probability FA,ri- By taking a £ {0.0, 0.1, 0.2, . . . , 0.8}, we obtain 7920 regulatory networks, 
20 for each triplet {p,a,r]) £ {0.2,0.4,0.6,0.8} x {0.0, 0.1, . . . , 0.8} x {0, 0.1, 0.2, . . . , 1}. 
Finally, for each one of those dynamical systems we run 30 initial conditions randomly 
taken from [0, 1]^'^'', and we let the system evolve in order to determine the corresponding 
oscillatory subnetwork and the asymptotic period. 

The characteristics of the oscillatory subnetwork depend strongly on the parameter o, and 
less noticeably on r]. We have observed that the mean size of the oscillatory subnetwork 
decreases with a in a seemingly monotonous way, and does not appear to depend on r/. 
The mean number of connected components never exceeded 2. A crude inspection of the 
degree distributions of the oscillatory subnetworks suggests that they form an ensemble 
of the Erdos-Renyi type, nevertheless, any conclusion in this direction requires a more 
detailed numerical study. Finally, the mean value of the asymptotic period seems to grow 
with both a and r]. It increases monotonously with a, while for each fixed a, it reaches a 
local maximum around r] = 0.5. 

As a matter of illustration, in Figure 1 we show two digraph with vertex set V = 
{1, 2, . . . , 50}. At the left side we draw a random digraph Q £ (dy according to the 
distribution ¥p, with p = 0.1. The signs of the arrows correspond to a random choice 
according to Pa,??, with r] = 0.5, i. e., both signs have the same probability to be chosen. 
The digraph at right is the oscillatory subnetwork Q = Gosc{Fg,i7,T,a,^), corresponding to 
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a certain choice of thresholds T € [0, 1]"^ and initial condition x G [0, 1]^. The contrac- 
tion rate was a = 0.2. In the oscillatory subnetwork, the direction of an arrow (n, v) is 
indicated by marking head v with a x sign. 




Figure 1. An Erdos-Renyi network and corresponding oscillatory sub- 
network. The underlying graph was built with a connection probability 
p = 0.1, and the oscillatory subnetwork corresponds to a random initial 
condition, with random threshold values, and contraction rate a = 0.2. 
Signs are codified by colors: blue for and red for -1. They appear with 
equal probability rj = 0.5. In the oscillatory subnetwork, the direction of 
an arrow is indicated by marking its head vertex with a x . 



4.3. Barabasi Albert. 

Following the procedure described in 4.1 we generate, starting from a complete graph 
in rriQ = 5 vertices, ensembles of random digraphs of the Barabasi-Albert kind of size 
N = 100. For each r/ G {0,0.1,0.2, . . . , 1}, we generate 20 random graphs, then we turn 
them into directed graphs by choosing one of the two possible directions with probability 
1/4, and both directions with probability 1/2. Finally we randomly assign thresholds, 
and signs with probability P^,^, to each one of the 20 resulting digraphs. By ranging 
a € {0.0,0.2,0.4,0.6}, we obtain 880 regulatory networks, 20 for each couple (a, ry) € 
{0.0,0.2,0.4,0.6} X {0,0.1,0.2,... ,1}. Then, for each one of those dynamical systems, 
we randomly select 30 initial conditions in [0, 1]^'^'', and let the system evolve in order to 
determine the corresponding oscillatory subnetwork and the asymptotic period. 

For this family of random networks, the mean size of the oscillatory subnetwork behaves 
in a similar way as for the Erdos-Renyi family, i. e., it increases nonmonotonously with r/, 
with a local maximum at r/ = 0.5, and decreases monotonously with a. The mean number 
of connected components (ranging from to 10) shows a similar behavior with respect to 
both T] and o, with a very noticeable increase around r/ = 0.5. On the other hand, the 
distribution of the asymptotic period becomes broader as we approach rj = 0.5, and for 
each value of ry, the mean value of the asymptotic period grows with a. It is worth to 
notice that these distributions become broader as a increases, and at the same time the 
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number of connected components grows. This behavior is consistent with the prohferation 
of periods predicted by Corohary 1, in the case of a statistical ensemble whose oscillatory 
subnetworks have many connected components. 

In Figure 2 we show two digraph with vertex set V = {1, 2, . . . , 50}. At the left side we 
draw a random digraph Q (z &v according to the distribution Pba- The sign of the arrows 
is randomly chosen by using the distribution Pa, 77 with rj = 0.5, i. e., both signs have the 
same probability. The digraph at right is the oscillatory subnetwork Q = GosciFg,cT,T,a,^), 
determined by a random choice of thresholds T € [0, 1]^ and initial condition x € [0, 1]^, 
with contraction rate a = 0.2. In the oscillatory subnetwork, the direction of the arrows 
is indicated by marking the head with a x . 




Figure 2. A Barabasi-Albert network and corresponding oscillatory sub- 
network. The underlying graph was built according to the procedure de- 
scribed in Paragraph 4.1, starting with niQ = 3 vertices. The oscillatory 
subnetwork corresponds to a random initial condition, with random thresh- 
old values, and contraction rate a = 0.2. Signs are codified by colors: blue 
for +1, and red for -1. They appear with equal probability r] = 0.5. In the 
oscillatory subnetwork, the direction of the arrows is indicated by marking 
the head with a x . 



In Figure 3 we show the behavior of the probability to approach a fixed point, and the 
probability for the asymptotic period to be equals 2, both as a function of a and t]. This 
picture summarizes the statistics of the 26400 orbits generated on Barabasi-Albert like 
networks, and presents the symmetry 77 1—r] predicted by Proposition 1. This symmetry 
is expected in an ensemble of regulatory network whose underlying digraphs have no odd 
length cycles. Though the underlying networks of the experiments we performed have 
both odd and even length cycles, there is a larger proportion of even lengths amongst the 
cycles of small length. 



5. Proofs 
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Figure 3. At the top, the probabihty to approach a fixed point in the 
Barabasi-Albert ensembles, as a function of a and rj. At the bottom, the 
probability for the asymptotic period to be equals 2, also as a function of 
a and r/, for the same ensembles. 

5.1. Proof of Theorem 1. 

Given Q = {V, A), we will use the distance dmax(x, y) := maxj |xj — yj| in both [0, 1]^ and 
[0, 1]^. As usual, dmax(x,S) := inf{(imax(x, y) : y e E} for x G [0, 1]^ and E C [0, 1]-^. 
First notice that for Q, a and a fixed, the set 

(10) Og,^,a := I (T, x) G [0, 1]^ X [0, 1]^ : inf d^^x (F* Jx), At) > o| 

is open. Indeed, if (T, x) G Og^a,a, then dmax (^^g a T a(^) ^ '^T^ ^ "^^i ^'-'^ some e > and 
all i ^ 0. The orbit does not change if we perturb T, i. e., for each T' G B^{T) we have 
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^g,a,T',a(^) = ^,cT,T,ai^) ^ G N. Furthermore, since dmax (^F^ ^^ j,, ^^{x) , At'^ > 2e, 

then dmax (^F^ ^^^t' ,aiy) ^ ^T'^ > ^ foi' all y € ^^(x), therefore ^^(r) x ^^(x) C Og^a,a- 
Since for each Q, a and a fixed, the set Og^o-.a is open in [0, 1]"^ x [0, 1]^, then 

(11) Oa,v:=(j \J U{(Fg,,,T,a,x) : (T,X) eOg,<,,4 

is measurable in Tay ^ [0) ll^- 
S.l.l. The asymptotic period is measurable in Oay- 

Fix Q^a and a, and take (T, x) G Og^a-,a such that P{Fg^^^x,a,^) = t- There exists 
y G Perr(Fg,CT,T,a) and iV € N, such that 



, ^ , initio C^max (i^g,^,T,a(x), At 
dmax (F^,,,T,a(x),y) < £ ■= 



Hence, for each i ^ we have dmax y^g^aTai^)-' aT aiy)j < same argument 

as in the previous paragraph, dmax (^Fg^!^, ^(y), Fg ^j,, „(y)^ < a^e + a^~^^£ < 2a*e for all 

t ^ 0. Therefore P{Fg^(j^T' ,a,y) = t for all (T",y) G BeiT) x -Be(x). In this way we have 
proved that for each Q, a and a fixed, the set 

41a ■■= {^x) G [0,1]^ X [0,1]^ : P(Fg,.,T,a,x) =r} 

IS an open set in [0, 1]^ X [0, 1]^, therefore 

P-\t} n Oay := U U U {(^e,.,T,a,x) : (r,x) G O^^i,} 

Q ^€{-1,1}-" 

is measurable in Tay x [0, l]^• 
5.1.2. The asymptotic period is finite in Oay- 

For (Fg^o-,T,a) x) G Oa,y let X G closure (^{-fg CTTa(^) • * ^ N}teN^ be such that 
inf dmax (^g,<T,T,a(x), At) = d^ax (x, At) = e > 0. 

Notice that (Fg;,^,T,a, x) G Oay too, and inft;>o dmax (-Pg,a,T,a(^)' At) = dmax (x. At) = e. 

For N such that d^^ax (i^^^^T.a (x) ' ^) < e we have d^^ax Jx), F*(x)) ^ a*e 

for all t ^ 0. Now, for r G N such that < 1/4 and dmax (^F^j^'^ ,j (x) , x) < e/4, we have 

dmax(FeV,T,a(x),x) ^ d^ax (^5';;+?, Jx) , x) + <ax (^5^;;+?,, (x) , (x) ; 

< e/4 + a^e<e/2. 
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Since inf^^o c^max (^-^g,cr,T,a(^)' ^Tj = e, then 

for all € N. Hence y := linifc^oo -^g ^ T a (^) exists and satisfies Fg ^Taiy) ~ y- Further- 
more, since dmax (x, y) < e Zlfcli 2-*^"^^ = 2e/3, then 

N+kr /\ ~\ ^ T f jpN+kr ' - -i-- '-^^ - '-^ ^^'^^ 



dmax (i^^:+^:.(x),yj ^ rfmax (F™ (x) , (x) j + d^ax yj < 

for each A: € N, which implies hmt_oo -Fg,cr,T,a(^) ~ ■^e,(T,T,a(y) = 0' where y = F^^^ j,^^{y) 
with n = — mod r. It follows from here that P{Fg^cr,a,T,^) = In addition, since 
y e closure ({-^g^^^T.aW ^ * ^ NltgNj, then mino:^t<r c^max [F^^^^^j,{y), At) = e. 



5.1.3. The complement of Oay has zero measure for a < maxugv'(ld(ii) + 1) ^. 

Fix g,a, a, and for (r,x) ^ Og^^^a let -D^*) := ^e,<7,a(i^g,^,T,a(^)) each t G N. Then, 
there exists a sequence {i^ € Nj^gN, and u,v €V, such that 

T(,,,) = ^lirn^ ^a*''x„ + (l-a)^a*-il)^-*)^ 
(12) = (1-a) lim lj2^'~'Dt 



For each x G [0, 1]^ fixed, define the fiber 

(13) Cg,„,a{^) := {T G [0, 1]^ : (r,x) ^ Og,^,a} . 

Then, according to Equation (12) we have 

Cg,a,a{^) C U |r G [0, 1]^ : r(„,„) G ^ X closure (J^g,a,«)| 

where 

^e,a,n := U E : (k.)l-o' e {0, 1, . . . , Id(n)}^ . 

Ar=l l s=0 J 



If a < (Id(u) + 1)~^ for some u ^ V, then closure (fig^a^u) is a Cantor set of dimension 
log(Id(?i) + l)/log(l/a) < 1. Hence, Fubini's Theorem implies, for a < max^gv'(Id(M) + 
1)~^, that the fiber Cg^o-,a(x) has zero ^^A-dimensional volume for each x G [0, 1]^. 
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5.1.4. The complement of Oay zero measure for a ^ max„gv'(Id(ii) + l)""*^. 

The fact that vol (Cg.o-,a(x)) = for each x € [0,1]"^ and arbitrary a € (0,1), derives 

from a result by Krughkov and Rypdal. Theorem 2 in [12] gives an upper bound for 
the topological entropy of the symbolic system associated to a non-degenerated piecewise 
affine map. The upper bound is related to the exponential rate of angular expansion under 
the action of the map. For piecewise affine conformal maps, which is the case of discrete- 
time regulatory networks, this upper bound vanishes. Hence, Kruglikov and Rypdal's 
theorem directly applies to our case, implying that 

(14) hm sup = 0) 

where for each Q,cr, and a fixed, 

4tIu) ■■= {(W.)f=o' € {0, 1, • • • ,Id(n)}^ : I)g,.,T(Fa,.,T,a(x))« = G 

Therefore, for each T, there exists N = N{T, x, n) G N such that #>cj^J^ ^ ((l+a)/2a)^. 
Fix u£V, and for each iV G N and P C {0, 1, . . . ,Id(n)}^, define 

7^ := |x = X]"*'^* • i'^{k+i)N~i,- • • ,kfc7v) G P Vfc ^ o| . 

With this we have, for the fiber Cgaai^) defined in (13), the inclusion 
(15) 

Cg,a,ai^) C U U U G [0, 1]^ : r(„,,) G ^ x closure(T^)| . 

( (N) 

Now, for each m G N, a prefix k G ( ^(t^xu) ) defines an interval 



t=o 



and we have 



closure(7x>) = U clos^'^^ (-^k) • 

m=l Vke©™ / 

Since vol(Ik) = vol(closure(/k)) = Id(ii)a'^^/(1 - a) for each k G (^\xxu)) ' ^^'^^ 



vol(closure(Tp)) ^ {#V)"'a 



rnmN Id(n) 



1-a 

for each m G N. Hence, vol (closure (T^))) = for each D C {0, 1, . . . ,ld(u)}^ such that 
^ ((1 + a)/2a)^ . Taking into account (15), Fubini's Theorem implies, that the fiber 
Cg,o-,a(x) has zero ^^-dimensional volume for each x G [0, 1]^. 
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Using once again Fubini's, we finally conclude, 

Pa,„(^a,V X [0, 1]^ \ Oay) ^ ^ ¥G{Q)^A,n<y I ^_ VOl (Cg,<,,a(x)) dx = 0. 

e;ee5v<76{~i,iH ^xe[o,i]^ 

□ 

Remark 1. In the statement of Theorem 1 we can replace the probability measure P^^^ 
by any other Borel measure P in J-ay x [0, 1]^ such that, for each Q € (J5y, a € { — 1, 1}"^; 
and X € [0, 1]^ fixed, the jj^ A- dimensional projection 

¥t{J) := P{(Fg,,,T,a,x) S J'ay X [0, 1]^ : T G J} 
is absolutely continuous with respect to the Lebesgue measure. 



Remark 2. From the arguments developed in Paragraph 5.1.3, it follows that for a < 
min^jgy (Id(ii) + the set Og^o-,a defined in (10) is open and dense. This ensures that 

orbits are generically uniformly far from the discontinuity set. As a consequence, periodic 
attractors not intersecting the discontinuity set are generic for a < min^jgy (Id(M) + 
On the other hand, the argument in Paragraph 5.1.4 ensures, as long as we consider 
distributions absolutely continuous with respect to Lebesgue for thresholds, that orbits are 
almost surely uniformly far from the discontinuity set. Hence, in the general case, almost 
all initial conditions converge to periodic attractors not intersecting the discontinuity set. 



5.2. Proof of Theorem 2. 

This theorem consist of two claims. The first one establishes the stability of the oscillatory 
subnetworks, while the second one gives the equivalence between the dynamics restricted to 
that subnetwork and the dynamics supported by the subnetwork considered as an isolated 
system. Both claims follow from direct computations. 

5.2.1. Stability of the oscillatory subnetworks. 

Since fa,ri{Qosc = G) > 0, then, there exists Q G ©y, a G {—1,1}"^, and (T, x) € Og^a,a, 
such that QosciFg = Q. As proved in Paragraph 5.1.2, there exists < e < 1/2, 

r G N, and y G Per^(Fg.^^^^^^), such that limt_,oo - F^^ _ ,^Jy) = and 

From this it readily follows that G osc[Fg, a, T, a iY) = Q for each (T, y) G B^{T) x -Be(y). 

Let us now extend Q, a, and T, and redefine y at each vertex v £ V such that Idosc(^') = 
Id(t;) (remind that Idosc(^^) := #{u G V : {u,v) G A = Aosc})- Choose u ^V, and include 
the arrow {u,v) in A. Then, if y^ < 1 — 2e, define cr(M^t,) = 1 and T(^u,v) = 1 — e. Otherwise, 
if Yu > 2e, make (T(„ ,;) = —1 and T(^u,v) = Finally, redefine y^, = ld{v) y^,/(Id(u) + 1) + 
l/(Id('y) + 1) at that vertex. A simple computation shows that, in the redefined system, 
OosciFg,a,T,a,y) = Q, for each (T, y) G Bs{T) x -Be(y). Therefore, taking into account that 
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by hypothesis fciG) > for all Q € ©y, we can assume without lost of generality that 
Idosc(f ) < Id(u) for ah v eV. 

By definition, if {u,v) ^ A \ A then ^(^^t,) := H {a f^u,v){y\i ~ ^(n,i)))) remains constant in 
time, as well as D{v) := J2{v u)eA\A^{v,u) for each v & V. Let us now modify a and 
T for the arrows leaving Q. For u G V, and each v ^ V, D{u) > 0, make Tf^^^^-^ i— s- 
T{u,v) X D{u)/l<l{u), and 

^"'"^ I -^^(n,^,) if 0{u,v) + (^{u,v) G {0, 1}. 

Otherwise, if D{u) = 0, then T(„ „) i-^ T'(„ x (Id(n) - Idosc(ii))/Id(n) + Idosc(^i)/Id(n), 
and 

-0"(n,i)) if (^{u,v) + 0"{n,D) G {-li 2} 
if ^{u,v) + (^{u,v) € {0, 1}. 

This modification in a and T uncouples the dynamics on V from that on so that now 

we can set = a for arbitrary a G {—1, l}"^- In this way we obtain a digraph extension 
Gext '■= Q ^ Q, and for each a € { — 1, 1}"^ an extension Ucxt := o" G { — 1) 1}"^ with ^extU = 
a. For any of this extensions, a direct computation shows that limt^oo ^ ^ t T S^"" ~ 
D{v)/ld{v) = y„ for v^V, and ,r Jx)„ G [D{v) /U[v),{D{v) +MZXv))mv)], 

for all t ^ and f G whenever 



'{u,v) 



X G J X J' := I ]J [y„ - e,y + e] n [0, 1] I x | J] 

and 



^ i^(^;)+Idosc(^) 
Id(w) ' Id(w) 

Z^(n) + Idosc(n) 



Id(u) ' Id(u) 



T G ZxT' := J] [r(„,„) - f + e] n [0, 1] X J] 

\{M,i.)eA\A / \(«,»;)gA 

Therefore, ^osc(i^g,,,,s,,,,r,a) x) C g for each (T, x) G T x J := (T x T') x ( J x J^'), and 
aG{-l,l}^. 

5.2.2. Equivalence to the subnetwork considered as an isolated system. 

We will use the same notation as in Subsubsection 5.2.1. Let a := ^exti and each (T, x) G 
(T X T) X (J X J') , and t G N, let x* := F| _ ^ (x). For each v €V, the term 

-D(t') := X] (^(«-'') ~ '^m)) = ^ ('^(«-^) (^g,<7,f ,a(y)« ~ '^M 

{u,v)eA\A {u,v)eA\A 
remains constant in time. Since = (T, 

M t , Idosc(^) / (1 - «) V- rrf- ft \\\^(^ n)RM. 



) ' 
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which can be rewritten as 

(16) $y(x^+^), = a^vi^% + ~ ° Yl ('^K.) («'y(x*)„ - «>a(T)(„,,))) , 

OSC \ / y \ _ 7 

(u,v)eA 

where $a : ^ x Z' ^ [0, 1]^ and $y : J x J' ^ [0, 1]^ are affine transformations defined 
as follows. We have ^aCT x T') = V^f + C^, with 

~\ Id(n) - J m \ -^('") 

^^^^)("'^) = id;::R^^("'^)' = "id;::(^' 

for each (u, v) S ^, and ^I'y (x x x') = VyX + Cy, with 

Id(u) _ , , Z)(i;) 

(P^xju = — — -- X x„, and [Cyjy - 



Idosc(^') Idosc(t^) 

for each v (zV. The result follows from Equation (16), which establishes 

Fg,a,'i>Anai^v{^)) = ^y(F^^^,,^^^,,T,a(x)) 
for all (T, X.) £ 2 X J', and from the fact that and $\/ are affine and surjective. 



□ 



Remark 3. In the statement of Theorem 2, instead of the hypothesis ¥g{G) > for all 
Q € (Sy, we could have directly supposed that the oscillatory subnetwork Q := {V^A) is 
such that Idosc(^^) < Id(t;) for each v . By assuming this, we avoid the condition V C^V 
as well. This alternative formulation allows us to consider the emergence of modularity in 
Barabdsi- Albert random networks, for which Pba(^) = for some digraphs Q £ <5v- 



6. Final Comments and Conclusions 



From our point of view, one of the most important theoretical problem in regulatory dy- 
namics concerns the relations between the structure of the underlying network and the 
possible dynamical behaviors the system generates. We have already established bounds 
relating the structure of the underlying network, and the growth of distinguishable or- 
bits [16]. Here, following [22], we have considered regulatory networks with interactions 
and initial conditions chosen at random at the beginning of the evolution. Within this 
approach, individual orbits become elements in a sample space, and the characteristics 
of the asymptotic behavior can be considered as orbit dependent random variables. The 
structure of the underlying network, random in this approach, is encoded in a probabil- 
ity distribution over a set of directed graphs. In two interesting cases we have explored, 
the asymptotic oscillations concentrates on a relatively small subnetwork, whose structure 
depends on both the proportion of inhibitory interactions, and the contraction rate. We 
have proved that the dynamics observed on this subnetwork is equivalent to the dynamics 
supported by the subnetwork considered as an isolated system. We interpret this as the 
emergence of modularity. Modularity allows us to predict asymptotic periods in regulatory 
networks admitting disconnected oscillatory subnetworks. 
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We want to remark Theorem 1, which has important technical imphcations. It ensures, 
as long as we consider distributions absolutely continuous with respect to Lebesgue for 
thresholds and initial conditions, that a random orbit converges to a periodic attractor. 
Furthermore, this periodic attractor does not intersect the discontinuity set. We have also 
proved, in the case of small contraction rates, that the periodic attractors are generic. 
This can be deduced from the argument in Paragraph 5.1.3. 

As mentioned above, in order to determine the distribution of the asymptotic period, 
and the characteristics of the oscillatory subnetwork, a detailed numerical study of the 
examples presented in Section 4 is required. Heuristic computations could guide those 
numerical studies, and it is our intention to proceed in this direction in a subsequent 
work. 
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